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I.  IBEGRY  OF  SURFACE  ARD  70UME  SOURCSS  OF  VOSTICm 
IH  ACOUSTIC  FIELDS 

A.  HKERODOCTIOR 

When  a sound  been  passes  through  a vlecoue  fluid,  voirtielty  is 
generated  as  a 8econd>order  effect  and  streaming  results.  An  expres- 
sion for  these  direct  cxirrent  streams  has  been  derived  by  Sckart.^ 

His  equation  (2;^a)  is  based  on  the  hypotheses  that  b » V/3  -f  ^ 
and  the  dynamic  shear  viscosity,  * pv,  are  both  constant  for  the 
fluid  density  variations  vhich  accompany  passage  of  a sound  wave, 
and  that  the  curl  of  the  firat-ordor  particle  velocity  is  zero 
X ^ » 0)»  Here,  P is  the  density,  v*  s |*»/p  is  the  kinematic 
bulk  vlscoelty,  and  v is  the  kinematic  shear  viscosity* 

The  assumption  of  the  density  independence  of  U axid  P*  bears 
investigatlcm.  It  is  certainly  true  tlmit  In  general  one  does  not 
expect  P or  P*  to  be  constant  in  a sound  field,  but  it  is  not  clear 
vhat  effect,  if  any,  their  variation  will  have  on  streaming.  The 
effect  of  variable  P is  left  for  a future  study;  but  a discussion 
of  density  dependent  P'  is  included  here,  and  it  will  be  shown  that 
the  vorticlty  equation  is  Independent  of  such  variation. 

The  second  assuaptlw,  i^amely,  7 x » 0,  is  an  important  one* 
It  is  probably  permissible  in  the  case  considered  by  Eekart  because 
the  departure  from  Irrotatlonal  behavior  is  small  if  the  beam  is  com- 
pletely "immersed”  in  the  medium*  However,  if  the  sound  beam  makes 
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cont&ct  vlth  a solid  surface  (say>  a solid  oibstacle  vlthin  the  sound 
field,  or  an  enclosing  tube),  then  the  particle  velocity  falls  rapidly 
to  zero  In  a transverse  direction  within  a thin  boundary  layer  meas- 
ured by  the  vlscou-:  wavelength.  Here  viscous  forces  becone  relatively 
iDQportant,  a velocity  potential  no  longer  exists,  and  the  assumption 
7x  « 0 Is  certainly  net  valid.  Actually,  this  nonzero  value  of 
the  curl  becomes  most  significant  in  such  cases  and  Is,  In  fact, 

p 

necessary  for  the  occurrence  of  the  streaming  discussed  by  Baylelgh, 
Schuster  and  Matz,^  and  Schlichtlng.^  It  will  be  sho^n  that  there 
are  three  sources  of  vortlcltyr  a term  Identical  to  the  Eckart 

soxnree  of  vortlclty;  the  Bayleigh  term,  which  appeeirs  when  the 
acoustic  field  is  rotational  (principally  at  solid  surfaces  where  the 
relative  quiescence  of  the  medium  Is  described  by  the  "boundary  layer"); 
and  S|j,,  a third  term,  new  to  discussions  of  streaming,  which  has  a 
nonzero  value  only  when  the  acoustic  field  Is  rotatlonatl. 

Finally,  a distinction  Is  made  between  "surface  sources"  of  vor- 
and  "volume  sources"  of  vortlclty  (S^),  and  solu- 
tions are  cited  for  common  situations  In  which  the  surface  sources 
d(Kslnate  the  stresmlng. 


tlclty  and 


, aURCES  OP  TORTICm 


In  the  formulation  of  the  Navler-Stokes  equations  given  by 
Leased  a typical  ooag>onent  eqvtatlon  Is 

* 3/6x  [xv-3  * ^ ^ (||  4 |sj] 


(1) 
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and  one  Is  confronted  immediately  vlth  the  decision  of  whether  or  not 
X = “ ~ ^ '^1^'  (Ji  are  to  he  considered  functions  of  position 
and  time  (and  thereby  functions  of  density).  Here  Xls  the  dllata- 
tlonal  viscosity,  Is  the  dynamic  shear  viscosity  In  poises,  and 
Is  called  the  dynamic  btUi:  viscosity.  P Is  the  fluid  density, 
u Is  the  vector  velocity  vlth  ccmiponents  u,  v,  v in  the  rectangular 
coordinate  frame  x,  y,  z..  D/Dt  Is  the  hydrodynamic  operator 
(d/dt  4-  u*7).  P is  the  mean  external  pressure  on  the  volume 
element.  Both  Eckart  and  Rayleigh  have  considered  that  P and  p"  are 
independent  of  density  at  this  point.  In  %he  following  derivation  of 
a more  general  eqiiation  which  predicts  streaming,  ii  will  be  con- 
sidered constant,  but  the  effect  of  variation  of  p ' vlth  density  will 
be  considered  by  the  expansion 

P'  * Pq  + (dp‘/dP)^Pl  “ (la) 

where  p^  is  the  bxilk  viscosity  at  equilibrium  density;  is  the 

first-order,  acoustic,  density;  (dp'/<iP)  represents  the  evaluation  at 

o 

equilibrium  density.  !Ehen  the  vector  force  equation  Is 

p^(3u/&t)  + P^U*7U  * -VP  + (V3P+P')  V7*U  - PVX  (7><u)  + (7P*)f7*u)  (2) 

The  notation  and  u^  will  be  used  for  the  acoustic  density  and 
particle  velocity,  respectively,  and  ^ for  the  second-order  particle 
velocity,  the  time  average  of  which  gives  the  stream  velocity.  Then 
the  first-order  continuity  equation, 
api/3t  + pQ^'u^  - 0 


(2a) 
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and  the  first-order  force  equation, 

pQciu^/at  » - + (4/3Pq+^o)^^.ui  • jIqVxCVjwi)  (3) 

are  the  saae  as  In  Eekart's  developoent.  Hovever,  t;he  second-order 
force  eq\iatlan  becomes  (assuming  the  pressure  Is  a function  of  den> 
slty  only): 

pQ^^/dt+pj^aui/at+PoCu]^ . V )uj^-.c§7p2-c<3(  dc/dp  )q7(p|  )+((%i  Vfip  )q(  ^ )^Pl 

-HoVxVxUg  - ^J^Vx7x^l^ 

With  the  assumption  ^ « constant,  » (d^/dP)^P^  0.  Further 

simplification  restilts  by  Introduction  of  P^(du^/dt)  from  £q.  (3)>  so 

that  vhen  the  curl  of  Eq«  (4)  Is  taken,  and  only  the  steady  state  Is 

considered,  there  results  the  general  vortlclty  equation: 

* + (-1/Pq)  (4/3+Wo^o)VplxV(^p3^/^t) 

— ' 

^ 

In  vhlch^2  ° ^ and  > 7 x Uj^.  (it  Is  noted  that  even  If  the 

medium  falls  to  satisfy  Hooke's  lav  so  that  (dc/dp)^  has  a nonzero 

value,  this  term  Is  multiplied  by  v(p|},  so  that  vhen  the  vortlclty 

Eq«  (^)  is  formed  the  term  vanishes.  This  fact  supports  Eekart's 

omission  of  this  effect  In  his  study  ot  the  Influence  of  viscosity  on 

vortlclty.)  Equation  (^)  Is  a general  equation  vhich  predicts 

steady- state  streaming  on  the  basis  of  p being  considered  constant 

but  vlth  no  assumption  as  to  the  constancy  of  p',  and  for  any 

acoustic  field,  irrotatloxmil  or  not. 

If  an  Irrotatl^ial  sound  field  is  assumed,  1L  ” ^ * 0 and  vhat 

R T 
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romalns  Is  Sckart's  £q.  (2S&)»  where  he  uses  the  definition  b « 4/3  '** 
It  is  noteworthy  that  the  possibility  of  variation  of  bulk  viscosity 
with  density,  which  is  considered  here,  does  not  affect  the  sources  of 
vorticlty*  It  is  clear  from  Eq.  (^)  that  the  case  0 is  a natural 
one  to  consider  when  one  is  Interested,  as  Eckart  was,  in  determining 
the  bulk  viscosity  . 

To  discuss  the  approach  of  other  workers  In  streaming,  it  is 
convenient  to  recast  Eq.  (^)  in  terms  of  the  kinematic  coefficients  v 
and  V*  » Observing  that,  for  constant  d''/dP  « 

Eq.  (^)  becomes 

= + (l/p^v^)(dy/dp)^  (4/3  + V */v  )^Pj^x^9p3^/9t 


+ (-1/vq)7x  + (lAo)(dv/»lp)o^  x(PiV3®l) 


(6) 


Rayleigh,^  and  Schuster  and  Natz^  started  with  Eq.  (2)  in  which 
as  well  as  ^ was  liiq;>licltly  considered  as  constant.  But  then,  in 
taking  the  curl  to  form  the  vorticlty  equation,  they  considered  a ro- 
tatlonal  acoustic  field  and  took  v and  V constant.  (Actually, 
Schuster  and  Matz  put  v ’ « 0. ) It  is  clear  that  only  the  term 
would  remain  after  these  assumptions,  and  this  is  precisely  the 
driving  term  which  Schuster  and  Matz  obtain  when  they  re>solve  the 
Rayleigh  problem  of  streaming  In  a tube,  using  vector  notation.  The 
assumptions  of  ^ constant  and  '*  ° ^/P  ■=  constant  are  equivalent  to 
assuiQ>lng  P e ccmstant  for  the  viscous  stress  terms,  only.  In  both 
works,  the  authors  apply  their  results  to  a compressible  fluid 


In  vhlch  P|f  0.  Schuster  and  Natz  raport  an  experimental  check 
vlthin  about  20  per  cent  of  theoretical  predictions  for  the  axial 
stream  velocity  halfway  between  loop  and  node  In  a standing  wave  In 
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an  air '■filled  cyl  Indr  Icsl  tube* 

2 •a 

In  addition  to  these  two  cases  of  streaming  In  a tube,  a 

search  of  the  literature  reveals  one  other  type  of  streaming  problem. 

Involving  oscillatory  motion  of  the  fluid  at  a solid  surface,  vhlch 

has  been  solved;  namely,  a rigid  cylinder  vibrating  In  static  fluid, ^ 

Sc^llchtlng  has  made  the  assumptions  v « constant,  v*  ^ o and  an 

Incompressible  fluid;  consequent  S and  6 would  be  zero,  and  the 

S T 

only  term  which  could  contribute  is  1^.  (Note:  Schllchtlng's  de- 

velopment which  leads  to  a streaming  velocity  does  not  parallel  the 
usual  approach  described  hers  and  in  ether  references  r He  obtains 
experimental  qualitative  agreement  with  his  theoretical  solution*) 

It  Is  noted  that  Sp,  the  third  term  of  Eq*  (^),  Is  new  to  dis- 
cussions of  strecualng  phenomena*  It  does  not  appear  when  one  assumes 
that  the  acoustic  field  Is  Irrotatlonal  or  If  one  assumes  that  P Is 
constant  In  the  viscous  stress  terms  during  passage  of  a sound  wave* 

In  no  case  In  the  past  have  the  authors  made  a critical  evalua- 
tion of  the  magnitudes  of  the  different  sources  of  vortlclty, 

although  It  is  clear  that  the  assumption  V x u^  0 Is  essential  to 

2-4 

all  three  solutions  In  which  the  acoustic  field  was  In  contact  with 
a solid  boundary.  IVo  additional  specific  cases  have  been  considered 
by  the  author;  the  streaming  caused  by  (A)  a plane  progressive  wave 

as 

in  which  ^ . « tube  radius  <X . and  (B)  a small 

viscous  acoustic 

rigid  fixed  spdxere  in  a plane  progressive  wave*  The  solution  to 
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problem  (A)  Is  presented  In  the  appendix.  At  this  tlme>  however^  It 

Is  possible  to  point  out  a distinction  vhlch  should  be  made  between 

surface  sources  of  vortlclty  and  volume  sources  of  vortlclty.  The 

descriptive  title  "surface  sources"  Is  used  for  all  sources  of  vor* 

tlclty  which  dependi  for  their  existence^  on  Vxu^^fOCi.e.^ 

and  Srp)|  and  the  title  "voltmie  source"  for  the  source  of  vortlclty 

which  exists  even  when  the  acoustic  field  Is  assumed  to  be  Irrotatlonal 

In  character  (l.e.  ^).  The  significance  of  the  terminology  Is 

this:  The  terms  ^ and  Sp,  depending  as  they  do  on  V x beccaae  ex> 

tremely  strong  In  the  neighborhood  of  solid  surfaces  where  the  viscous 

forces  are  relatively  la^rtant,  a potential -type  field  does  not 

exist,  and  large  valuec  of  7 x u^  are  generated.  On  the  other  hand, 

Sg  and  are  probably  negligible  in  regioxis  many  viscous  wavelengths 

distant  from  the  surface,  where  the  approximation  of  an  Irrotatlonal 

-+ 

acoustic  field  (l.e.  7 z u^  » 0}  is  generally  valid.  !Bie  term  Sg, 
however,  will  have  a zxonzero  value  In  the  medium  away  trott  the  surface. 
Further  support  for  this  terminology  lies  In  the  observation  that  Sg 
la  a function  of  tne  bulk  viscosity;  ^ and  are  not. 

Solutions  to  problems  (A)  axid  (B)  mentioned  above  clearly  Indi- 
cate that  In  these  two  cases  the  surface  sources  of  vortlclty  dominate 
the  streaming.  The  flow  velocity  caused  by  th<;  surface  sources  of 
vortlclty  Is  many  orders  of  magnitude  greater  than  that  caused  by  the 
voliB&e  soxirces  for  these  two  cases.  Althou^  it  has  been  suggested 
by  Eckart^  that  the  source  7ppX7dpp/dt  mlg^t  explain  the  stream- 
ing at  the  mouth  of  the  Helmholtz  resonator,  the  resiiilts  of  Schuster 
and  Matz  for  a standing  wave,  and  the  author  *^s  calculations  for  a 
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i 

I . T 

i 

traveling  vave  vhlch  fills  a tube  seem  to  liqply  otherwise.  Possibly  j 

, i 

the  theoretical  solution  to  the  problem  of  streaming  at  an  orifice  ! 

In  a plate  vhlch  blocks  a standing  vave  in  a tube^  Is  also  determined  | . 

i i 

by  the  surface  sources  of  vortlcity.  In  general^  the  final  assign-  J 

ment  of  the  Importance  of  the  source  terms  can  be  made  only  after  an 
csalxiation  of  the  terms  of  Eq.  (^)  through  a knowledge  of  the  actual 

first-order  acoustic  field.  The  difficulty  of  solving  Eq.  (^)  for  , | 

i 

i 

any  but  the  simplest  acoustic  field  and  the  simplest  geometry  Is  ^ 

obvious.  j ' 


! 

i 
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II.  SEBEAMma  CAUSES}  B7  A SOUHD  BEAM  IN  A CLOSED  TDBE 


The  analyslB  of  the  preceding  section  provides  a convenient 
hreedcdovn  of  vortex  streaioing  phenooena  in  tvibes  ac  due  to  surface 
sources  of  vorticity  which  predominate  in  a guided  acoustic  wave  or 
volis&e  sovtrces  which  are  most  important  in  an  acoustic  beam.  Ob- 
viously, the  latter  problem  is  cf  more  interest  since,  as  pointed  out 
by  Eckart,  it  appears  to  offer  a direct  approach  to  the  determination 
of  bulk  viscosities.  It  seemed  to  the  author  that,  preliminary  to  the 
use  of  the  streaming  eqmtlons  developed  by  Eckart  for  this  purpose, 
a complete  survey  should  be  made  to  study: 

a)  the  variation  of  streaming  with  the  pressure  amplitude  and 
distribution  of  the  acoustic  field,  and 

b)  the  variation  of  streaming  with  the  medium. 

The  dependence  of  streaming  on  distribution  of  acoustic  pressure 
could  be  studied  In  any  medlvos.  Then  the  use  of  a monatomic  gas  such 
as  argon,  with  known  values  of  bulk  viscosity  (zero),  shear  viscosity, 
and  heat  conductivity,  serves  to  "ccJ.ibrate"  the  experiment,  so  that 
the  constant  "b"  may  be  determined  and  interpreted  for  other  media. 

To  accon^llsh  these  goals,  it  is  iMcessary  to  knew  the  absolute 
values  of  the  acoustic  pressures  and  the  absolute  values  of  the  stream- 
ing velocities  at  all  points  over  a cross-section  of  a tube  filled 
with  a gas  of  known  purity.  The  equipment  used  is  described  in  section 
II-A.  The  extension  of  streaming  theory  based  on  the  observed  field 


and  comparison  of  predicted  streaming  Telocltlee  vlth  observation  Is 
taken  up  In  section  II>B.  Determination  of  bxilk  viscosities  of  moist 
air  and  dry  nitrogen  after  calibration  vlth  argon  Is  discussed  In 
section  II-C. 


A.  EQUIFMEOn!  AND  EXPSHMSfiTTAL  EROCEDUBES 

The  ultrasonic  streaming  tube  used  vas  a glass  pipe  of  six  Inches 
diameter  and  four  foot  length.  The  movable  source  holder  at  one  end 
of  the  tuba  held  a barltop  titanete  disk  transducer.  The  tube  vas 
terminated  by  a mcvable,  glass-vool  covered  end  piece.  Sxe  over- 
all appearance  Is  shown  In  Figure  1. 

The  source  vas  a one-half  Inch  thick,  three  Inch  dleuaeter  disk 
of  pre-polarlzed  barium  tltanate  ceramic.  It  vas  supported  at  three 
points  In  its  oedlan  plane  and  vas  capable  of  being  driven  In  any  of 
several  modes  betveen  the  frequencies  1^0  EC  and  210  ED.  Since  the 
theoretical  solution  Is  based  on  the  eissumptlon  of  a radially  sym- 
metric field,  this  vas  checked  first  by  sprinkling  lycopodium  powder 
over  the  disk  during  vibration.  The  pattern  observed  for  the  mode 
selected  for  this  experiment  Is  shown  at  the  top  of  Figure  2.  This 
motion,  vhlch  Is  not  a piston-type  mode,  corresponds  to  a frequency  of 
185  KC. 

Having  determined  that  the  source  pattern  of  vibration  was 
radially  sysnetrlc.  It  vas  then  necessary  to  make  the  crystal  axis 
collnear  vlth  the  tube  axis.  An  optical  method  vas  used.  A small 
plane-parallel-faced  mirror  vas  pressed  against  the  radiating  face  of 
the  barium  tltanate.  Then,  olgliting  through  a set  of  Illuminated 


ACOUSTIC  STREAMING  TUBE 


LEFT  : Ba  Ti  O3  source 

RIGHT:  Glass  wool  absorber 

CENTER:  A.D.R  Probe  microphone  and 
preamp,  with  recorder 
and  drive 

FIG.  I 


l-ycopodium  pattsrn  of  sourco 

BOTTOM;  Probe  microphone  with  A.D.P, 
crystal  at  lower  left 


FIG.  2 
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cross-hairs  centered  on  the  axis  and  at  the  other  end  of  the  tube, 
the  crystal  vas  adjusted  until  the  original  cross-hairs  and  their 
iaage  in  the  mirror  vere  superimposed  at  the  center  of  the  mirror. 
Estimates  of  errors  lead  to  the  concluslcn  that  this  sii'^llfled 
Oausslan  eyepiece  technique  would  give  a angle  of  dewlatlon 

of  1^  minutes  between  the  perpendicular  to  the  crystal  face  and  the 
tribe  axis.  This  amounts  to  the  center  of  the  beam  being  off-axis 
by  a maximum  of  slightly  over  one  millimeter  for  the  conditions  used 
In  most  of  the  experiment*  That  the  radially- symmetric  source  vibra- 
tions gave  a radially- symmetric  acoustic  field  when  aligned  as  de- 
scribed, Is  shown  by  the  field  recordings  In  Figure  3*  method  of 
obtaining  these  records  Is  described  below. 

Two,  l/8th  Inch  O.D.  probe  microphones  (lower  object  of  Figure 
2}  were  used,  one  movutted  transverse  to  the  tube  on  twin  guide  rails 
with  rack  and  pinion  drive,  the  other  mounted  through  the  end  plate, 
with  prefbe  axis  parallel  to  the  tv^  axis*  The  transverse  microphone 
was  driven  by  the  motor  of  a Bruel  and  EJaer  Recorder  and  was  used  to 
probe  the  acoustic  field  along  a diameter.  Different  diameters  could 
be  probed  by  rotating  the  soxirce  holder  about  Its  own  axis*  Records 
for  four  of  these  diameters  at  a distance  of  ^0  cm*  from  the  source 
are  superimposed  in  Figure  3*  These  records  show  the  essentlelly 
radially- symmetric  character  of  the  field*  The  end  microphone  was 
used  to  monitor  the  acoustic  field  during  streaming  velocity  observa- 
tions when  the  transverse  odoropbone  was  withdrawn*  The  sensitive 
element  for  each  microphone  was  an  ammonium  dlhydrogen  phosphate 
(A*D*P.)  longitudinal,  expander  beur,  approximately  0.8  mm*  x 2 mm* 
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X 7 (see  bottom* left  corner  of  Figure  2)  held  at  Its  midpoint 
so  as  to  protrude  sll^tly  from  the  isrohe  tube.  A double  shielded 
cable  connected  the  crystal  to  a preaiqpllfler  adaptor..  The  sensitiv- 
ity of  the  microphones  vas  determined  by  direct  coaparlson  vlth  a 
W.E.  6^0AA  condenser  microphone  which  vas  calibrated  by  reciproc5.ty 
to  an  accuracy  of  + 0.3  db. 

Thermal  convection  proved  to  be  a most  troublesome  Interference. 
To  minimize  such  deleterious  effects  it  vas  necessary  to  place  the 
streaming  tube  in  the  small  U.C.L.A.  anecholc  chamber,  which  is  in 
turn  within  a temperat\ire  ccmtrolled  1°C.)  ro<xn.  Since  body  heat, 
itself,  caused  thermal  Btreamlng,  motion  of  tobacco  smoke  or  Lyco- 
podium tracing  particles  vas  observed  by  telescope  throu^  a hole  in 
the  wall  of  the  anecholc  chamber.  Line  illumination  of  the  particles 
In  the  tube  vas  provided  by  a line  filament  lamp  with  spherical  lens 
vhlch  gave  a line  image  at  a distance  of  about  one  foot  fr«n  the 
ll^t  sovorce.  An  infra-red  filter  was  effective  if  the  ll^t  vas  net 
used  for  more  than  30  minutes  at  a time.  Heating  of  the  pre-amplifior 
and  the  probe  microphone  placed  a 1^  minute  operation  limit  on  the 
microphone.  Heating  of  the  source  transducer  Inposed  the  most  severe 
restriction  on  the  experiment.  After  only  three  minutes  the  source 
would  be  warm  enough  to  drive  thermal  convection  streams  for  hours. 

The  telescc^  vas  provided  vith  a reticle  vhlch  vas  then  cali- 
brated for  fields  of  view  at  various  distances.  In  operation,  the 
position  of  the  lllumli^ted  line  was  set  by  coving  the  light  source 
holder  by  means  of  its  support  rod  vhlch  extended  outside  the  anecholc 
chaaiber.  The  telescope  vas  then  focused  on  the  particles  which  ve^’a 
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llluBlnated  by  the  li^t|  and  a stop  clock  vhlch  read  to  hundredths 
of  seconds  was  used  to  time  the  passage  of  the  particles  past  a cali- 
brated reticle  distance* 

The  electronle  circuitry  Is  shoim  in  the  block  diagram  of  Fig- 
ure The  barium  tltanate  transducer  was  driven  by  a signal  from  a 
GJR,  72£a  oscillator,  vhoso  output  vas  aapllfled  by  the  amplifier 
section  of  a W.E.  17B  oscillator.  The  microphones  were  adapted  to 
the  Vestern  Electro-Acoustic  Pre-angpllf  ler  and  Microphone  Complement 
100b,  with  bias  voltages  blocked  from  the  A.D.P.  crystals.  The  out- 
put of  the  C<miplement  vas  sent  through  a Navy  B/iK  7 receiver  for 
amplification  and  filtering.  This  output  vas  applied  to  the  volt- 
meter or  to  the  Bruel  and  EJaer  recorder. 

The  rather  cotaplex  source  holder  shown  in  Figure  1 vas  neces- 
sary to  provide  for  alr-tlght  movement  in  the  (unfortunately)  slii^tly 
tapered  glass  tube*  Since  all  experiments  vere  carried  on  at  atmos- 
pheric pressure,  only  diffusion  at  no  difference  of  pressure  could 
bring  In  lng)urltles.  To  minimize  such  leaks,  "0"  rings  vere  used 
throughout.  The  source  holder  slides  on  two  '*0'*  rings  vhich  are  kept 
tight  by  twelve  bolts  passing  throu^  their  own  "0”  rings.  Tightening 
the  bolts  squeezed  the  large  ”0''  rings  out  to  contact  the  variable 
I.D.  tube.  The  microphone  probes  also  passed  through  ”0"  rings.  One- 
quarter  Inch  steel  end  plates  supported  the  pound  force  of  air 
pressure  exerted  on  the  source  holder  and  end  piece  during  evacuation 
for  change  of  gas.  The  equivalent  leakage  area  vas  calculated  to  be 
10'2  mm^  by  the  change  of  pressure  as  air  seeped  Into  the  tube  after 
evacuation.  Diffusion  throuc^  sxich  an  area  vould  give  only  about 
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0.01^  laiipurlty  In  100.  hours.  No  gas  sao^le  was  ever  kept  In  the 
tube  for  more  than  24  hours  In  the  experiment. 

The  gases  used  In  the  experiment  vere  moist  laboratory  air, 
cGdoserclal  argon,  and  ccsunercled,  oil-pumped  nitrogen.  The  nitrogen 
was  asserted,  by  the  sxqppller,  to  be  99*7^  pure  with  0.01^  water  mole- 
cules. Both  Lycopodium  powder  and  tobacco  smoke  vere  used  as  tracing 
particles,  with  no  real  difference  in  the  streaming  observed.  It  is 
felt  that  the  water  molecule  Impurity  was  no  more  than  0.03^  and  had 
no  measurable  effect  on  the  streaming  in  nitrogen  or  argon. 

B.  mPESmSCE  0?  £nS£AMING  ON  ACOUSTIC  FIELD 


Using  only  the  first  term  of  equation  (^),  Eckart  shows  that  a 
radially- syissstrlc  progressiva  sound  beam  of  acoustic  pressure, 
p^  * ^1°!  * ■ wt) 

In  a tube  of  radius  r^,  will  cause  streams  which  have  a distribution 
with  r given  by: 

Ug(r)  = r (8,r)P^(s)ds  + P(r^^  - r^)  (?) 


where  K » (^  bk^)/(P^c^) 

k « 2^tA 

w B 2^f 


r 


P 


B static  fluid  density 
*=  speed  of  sound 


V3  + 


» log  (^yr) 


When  8 


(8) 


s-lc 


-T  r\9  ( /«  ^ 


when  s - r 


4.  ,**0  , 2 3.  _2. 


- (K/r_  )/^  (rro  -r  ) P (r)dr 


(9) 
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The  applicability  of  these  equations  to  a real  acoustic  field 
In  air  has  been  experlnentally  studied  by  the  author  vlth  regard  to; 

1)  varlatlcm  of.  streaming  vlth  acoustic  Intensity, 

2)  variation  of  streaming  (both  on  axis  and  off)  vlth  acoustic  field 
distribution,  and 

3)  attenuation  of  axial  stream  velocity  near  a source  of  sound. 

1)  Variation  of  Streaming  vlth  Acoustic  Field  Intensity 

Equations  7-9  predict  that  the  streaming  velocity  should  be 
proportional  to  the  square  of  the  acoustic  pressure  at  a point  (for 
a fixed  field  distribution).  It  Is  necessary  to  investigate  this 
functioned  dependence  for  the  range  of  Intensities  to  be  studied  be- 
fore using  the  theory. 

Observation  of  tobacco  smoke  or  lycopodium  povder  tracing 

particles  vas  used  throughout  the  experiment  and  gave  accurate  values 

up  to  speeds  of  about  1.3  cm/sec.  Within  the  range  up  tc  1.3  cm/sec,  u 

2 

Is  apparently  proportional  to  p^  . This  Is  shcvn  In  Figure 

This  result  may  be  compared  vlth  tvo  criteria  vhlch  have  ap- 
peared In  the  literature.  First,  for  direct  current  flov  In  a smooth 
pipe,  the  critical  Beynold's  Number,  * v>2<lA/  vhlch  turbulence 
sets  In  Is  2000  or  greater  depending  on  entry  conditions  of  the 
stream.*^  Such  a criterion  applied  to  this  experiment  vould  predict 
onset  of  non-linear Ity  at  a flov  velocity  of  200  cm/sec.  Cn  the  other 
hand,  Llebermanu^  has  conducted  em  experiment  similar  to  this  one  (but 
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in  liquids)  and  found  that  the  constant  of  proportionality  between  Ug 

p 

and  p^  changed  at  a Reynold's  Nunber  of  about  20.  If  this  criterion 
held  here,  non-linearity  would  be  expected  at  velocities  greater  than 
U2  2 ^ ca/sec.  No  such  change  of  linearity  was  observed  in  the  range 
of  axial  velocities  from  zero  to  1.3  cm/sec. 


2)  Variation  of  Stre*nn<ng  with  Acoustic  Field  Distribution 

If,  for  simplicity,  one  assumes  an  *^ideal''  acoustic  beam  given  by 


Pl(r) 

\(r) 
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r - r. 


(10) 


where  r^  is  the  radius  of  the  non-dlvergent  sound  beam,  there  results 
from  equations  (7)-(l0)  a prediction  of  an  axial  stream  velocity,  for 
such  an  "ideal"  field: 

ug(0)  * 


r 2/  2 , 

, 2 2 ) '1  /'o  - 1 . 

I KPj^OTi  I 2 - log  ^i/^o 


SPo'c** 


(11) 


where  G « {r2^/r^^-l)/2  - log  r^/rQ 
b * k/3  ♦ v'/v 
I - Pi^(2PoC„) 

This  last  equation  was  derived  by  Eckart  and  used  by  Llebermann. 
Llebermann's  velocities  were  obtained  by  measurements  of  force  on  a 
diaphragm  which  had  been  calibrated  by  observation  of  tracing  particles 
moving  with  the  fluid. 

When  a probe  mlcro|dione  is  used  to  determine  the  actual  distri- 
bution of  acoustic  pressure  in  the  field,  it  is  possible  to  substitute 
this  Infonuitlon  into  equations  (7)*(9)  thereby,  to  get  a more 
realistic  prediction  of  the  streaming.  Without  question,  the  rctual 
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field  P^(r)  shoTild  be  used  to  calcxilate  the  streaolxig,  and  vbis  Is 
vliat  has  been  done  in  these  experlnents*  EoweTer,  It  Is  of  interest 
to  coaqpeare  the  axial  streaming  predicted  by  the  actual  field  with  the 
streaming  predicted  when  one  uses  the  "ideal"  field  assumption.  It 
will  be  proved  that  a gross  error  vould  be  made  in  the  streaming  pre- 
diction if  one  replaced  the  actual  field  used  in  this  experiment  by 
an  equivalent -energy^  "ideals"  field  as  described  by  equation  (lO). 

a.  Axial  Velocity 

For  the  special  position  r » 0,  but  with  the  general  field  dis- 
tribvrtloni  P(s), 

s - r and  (B,r)  * s log  r^/s. 

Then, 

UgCo)  » log  (^o/s)  i^(8)  - K/r^V  ^ (sri  - s^)P^(s)ds 

^ o o 

» q P^(r)  dr  (12) 

3 2 

where  Q » r log(rQ/r)  - r ♦ r /r^ 

Q(r)  is  a "veie^tlng"  feu:tor  which  determines  the  contribution  to  the 

2 

axial  velocity,  of  the  quantity  P (r).  Figure  6 is  a plot  of  Q(r)/rQ 
vs.  r/r^  and  emj^slzes  the  fact  that  it  is  the  distribution  of  P^(r) 
which  is  significant  in  determining  both  the  streaming  speed  and  di- 
rection on  the  axis.  (Note:  A tubular  field  with  a "hole"  on  the 
axis,  vould  result  in  axial  streaxaing  opposite  to  the  direction  of 
propagation  of  the  sound  beam.) 

In  the  following  discussion,  P^^  will  be  used  to  designate  the 
axial  (r  > 0)  maximum  value  of  -Uie  acoustic  pressure.  The  measure- 
ments at  72  cm  frcsn  the  somrce  are  considered  first,  and  the  axled 
maximum  acoustic  pressure  at  this  point  is  designated  by  ^1(^8 
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7 ehcftfs  the  average  (or  all  eight  radial  probings)  pressure  P(r), 

P^(r)  and  rP^(r). 

Ihe  stream  velocity  at  the  axis,  is  obtained  by  a gra^ical 

Integration  of  QF*(r)dr,  performed  in  Figure  6.  This  evaluation  gives: 

• ^o  ^107^  *0^ 


It  remains  to  answer  the  question,  "What  axial  streaming  is  pre* 
dieted  if  one  assumes  an  'ideal  field*  of  maximum  pressure  sj  1^,  where 
is  defined  so  that  the  force  of  its  radiation  pressure  against 
a hypothetical  end  diaphragm  eqimils  that  caused  by  the  actual  field?” 
For  both  fields  to  exert  the  same  fenree  ou  a diaphragm: 


2 \( 


\ 


1 \ ' 


^ /..  » I s.. 

r,vr/i  r or 


^0 

2 2 'o 

P ' 1 o 

P^(r) 

rdr-8^o/  r 
o 

Pl(r) 

dr 


where  r = "ideal”  beam  radius  ® r /2 
__  1 o 

As  defined,  this  mean  square  pressure,  , is  assumed  to  have  the 
above,  constant,  value  over  a "meas'urement  circle'"  equal  to  the  source 
area,  and  the  radiation  pressure  is  zero  everywhere  else. 

The  graphical  integration  is  performed  in  Figure  7 and  yields 
r 


^ o dr  - 0.0268  r^^ 


? tr\t  2\/  2^  2\  I 2 

C = \o/r^  ^10  ^ ^107 
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FinftXICTj  * 0«0085  (l^) 

Con^peirliag  equations  (13)  and  (l4).  It  is  clear  that  the  axial 
streaialsg  reloclt^  predicted  by  t^e  actual  acouatlc  field  Is  2,12  tines 
as  great  as  the  value  obtained  on  the  basis  of  an  equivalent  "Ideal” 
field.  Otherwise  stated.  If  the  axial  velocity  was  observed  and  then 
used  In  equation  (11)  to  solve  for  the  value  of  b,  the  value  of  b so 
obtained  would  be  approxltaately  twice  the  correct  value, 

b.  Off-axis  Velocities 

Not  only  the  magnitudes,  but  also  the  streaming  profile  (the 
variation  of  streaming  with  distance  from  axis).  Is  altezod  %/ht^u  the 
actual  field  rather  than  the  Ideal  field  Is  used  In  the  theory.  The 
method  used  to  get  an  analytic  expression  for  the  streaming  velocity, 
iigCr))  to  approximate  the  acoustic  field  P(r)  by  three  straight 
line  regions  defined  by: 

Beglon  A,  P(r)  « Pj^Car/r^  + b)  for  0<r  < rg 

Reglcn  B,  P(r)  ■ ^iQicr/r^  + d)  for  rgT  r<r^  (15) 

Region  C,  P(r)  ■ 0 for  r > r^ 

For  such  a field,  equations  (7),  (6)  and  (9)  provide  the  analytic 

e^gpxesslons: 

v^(r)  * {log  y ( y^/4(c2-a^)  + 2/3y3(cd-ab)  + j^/2(d^-b^)  ] 

o<r5rg  -log  z(c2*V*^+2/3  + a^/l6(y^-x^)  + 2/9ab(y^-x3) 

+ b^A  (y^-:^)  + c^/l6(z**-y^)  + 2/9cd(z^-y3)  -f  d^/4(z^-y2)  } 

+ prQ®(l-x^) 

ug(r)  « K PjQ^  ^®8  x[yVMa^-c^)  + 2/3y^(ab-cd)  + y^/2(b^-d^)  ] 

-log  z[Z  c^+2/3z3ca-*-d2/2  z^  + c2/l6(z^-x^)  + 2/9cd(z3-i3) 
+d2/4(z2-x2)}  + prp^(l-x^) 

ug(r)  - KPjjQ^r^^log  x[  y^/4(a2-c2)  + 2/3y3(ab-cd)  + y^/2(b^-d^)  + c2/4z^ 
r>^  + 2Acdy3  + d2/2z2)  + p (r/-r®)  (l6) 

P « - KPj^^  { -z^c2/6  + y^/6(-a^-Ks^)  - 2/5z^cd  ♦ 2/5y^(cd-ab)  + z^/k{c^-d^) 
+ y**/4(a®-b2-c^+d^)  + 2/3z^cd  + 2/3y^(ab-cd)  + d^z^/2  + y2/2(b^-d^)> 
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vhere  x ■ r/r^,  y ■ r^/r©  ajid  z - 

nie  averase  obsexved  pressure  distributions,  as  veil  as  the 
straight  line  approxlaatlons,  are  shown  In  Figure  9 for  three  dis> 
tances  from  the  source.  From  these,  the  following  numerical  values 
were  taken  for  the  two  oases,  26  cm  from  the  source  and  72  cm  from 
the  sovirce: 


constant 

Distance 

from  Source 

72  cm 

26  cm 

a 

-1,820 

-3.480 

b 

1.000 

1.000 

0 

-0.770 

-0.372 

d 

0.580 

0,223 

y 

O.UOO 

0.250 

z 

0.750 

0.600 

Use  of  these  constants  In  e<]uatlon  (l5)  then  gives  values  of 
U2(r)  from  which  the  curves  of  Figure  10  are  drawn  for  the  streaming 
profiles  at  26  cm  and  72  cm  from  the  soxiree.  The  streaming  pzadlcted 
by  the  "Ideal**  field  approximation  (for  beam  radius  rj^  • 
also  shown  In  Figure  10. 

All  three  curves  are  plotted  with  the  same  axial  velocity  to 
aid  In  coopering  the  profiles.  It  must  be  remembered  that  the  actual 
fle^d  method  would  give  streaming  magnitudes  different  than  those  ob- 
tained from  the  equivalent  "ideal"  field  approximation.  More  will  be 
said  about  this  in  the  next  section.  As  a chock  on  the  ability  of 
equations  (l^)  to  take  the  place  of  the  actual  field,  the  axial 
velocity  (for  72  cm  from  the  source)  obtained  here  Is  0.018U 
as  compared  to  0.0180  f obtained  by  iise  of  formula  Cl3). 

Figure  10  shows  that  the  stream  velocity  may  be  ejgpected  to  drop  off 


( K ^0^' 
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from  the  axial  value  more  rmpidly,  that  the  ohaafie  In  direction  ooctun 
nearer  to  the  axis,  and  that  the  eounter*flov  Is  relatively  slover, 
than  for  the  "Ideal*  field.  Iheae  dlfferenoea  are  greater  at  26  cm 
than  at  72  cm  from  the  source. 

fhe  plotted  points  In  Figure  10  are  axperlaantal  values  of  the 
streaalag  at  26  om  from  the  souroe.  Ihe  spread  In  the  values  at 
t/Vq  • 0.133  nay  he  explained  hy  an  en^r  of  only  2 lam  In  positioning 
of  the  light  sousroe.  It  may  be  said  that  ehen  the  observed  acoustic 
field  is  used  to  predict  the  stressing  velocity  profile  a ve:y  good 
cheok  vlth  experlmmt  Is  obtained. 

3 ) Quantitative  Yerlf Icatlon  of  Stream  Telocity  Prediction 

It  remains  to  check  the  namltude  of  the  streaming  velocity  ty 
study  of  a gas,  such  as  argon,  of  known  physical  constants.  inter* 
pretatlon  of  the  cause  of  acoustic  streaming  Is  tied  up  In  the  constant 
b of  eq^iatlon  (ll).  As  used  by  Bokart,  b m U/3  -f  » but  Sekart 

did  not  consider  thermal  ccnduotlvlty  or  relaxation  processes,  as 

9 

such.  More  recently,  Maxkhaa  has  extended  the  streaming  theory  to 
consider  the  effect  of  relaxation  processes. 

9 

Instead  of  using  the  realtlon,  pj^  ■ , one  takes  account  of 

the  possibility  of  acoustic  pressure  depending  on  the  rate  of  change 
of  density  by  expressing 

Pi  • (17) 

where  the  dot  Indicates  differentiation  with  respect  to  tJ.ne,  and 
where  is  some  kind  of  relaxation  constant. 
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lha  first-order  aquation  (3)  becoMs 
PoT?^  - - Co^Vpi  +bjio77.ui  7 X Pi 

O V 

But  since  the  aquation  of  continuity  (2a)  is 

®Pl^t  ♦ pQ  V'Uq  » 0 (2a) 

it  is  possible  to  vrita 

Po^^lM  - - Cq^Pi  pQ  77^1  + biJ^qy-Ui  x Ri  ■ - C^^Pi 

+ bVo’’*  - 1*0  7 X ?!  (18) 

vhare  b*  » b + (^P^/Pq  • V3  + »*4i^  +(Rpo^o 

The  argument,  vliioh  was  first  proposed  in  the  litexature  by  Fox 

10  Q 

and  Hertfeld,  and  has  been  strengthened  by  the  vork  of  Mixkham'  and 

V«  Hyborg,  ^ is  that  the  momentum  of  the  streams  is  evidence  of  the 
conversion  of  acoustic  momentum  lost  by  the  attenuated  sound  beam. 

According  to  this  eogument,  the  constant  b*  in  equation  (l6)  incor- 
porates all  attenuative  meohanisns.  Actually,  it  has  long  been  known 
that  even  shear  viscosity  and  thermal  conductlvi^  may  be  considered 
as  relaxaticn  phenomena,  the  relaxation  times  being  close  to  the  time 
between  molecular  collisions.  From  this  point  of  view,  thereal  con- 
ductivi^  may  be  included  in  b*  by  writing 

(r-1)  ^/Cp 

idiere  y ■ ratio  of  specific  heats 
X ■ thermal  conductivity 
Cp  ■ specific  heat  at  constant  pressure. 
'The  logical  conclusion  of  such  an  argument  is  that  one  may  write 
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D*  - - V3  ♦ tir-i)  «/CvCp)  + n'4i  C19) 

l£^H 

^ere  a ■ anplltude  attenu&tlon  coefficient, 
nepers  per  unit  length 
Ic  * 2ttA  * aootistlc  vave  number, 

and  that  the  determination  of  t*  either  through  a streaming  experiment 
or  through  attenuation  measurements  vould  give  the  same  results. 

In  equation  (19),  It  Is  clear  that  streaming  Is  proportional  to 
the  attenuation  and  Inversely  proportional  to  the  shear  viscosity  of 
the  medium.  That  Is,  the  streams  are  driven  by  the  rate  of  loss  of 
momentum  from  the  acoustic  uave,  and  are  retarded  by  the  shear  viscos- 
ity (at  the  vails  of  the  tube).  The  veil -known  classical  mechanisms 
of  shear  viscosity  and  thermal  conductivity  have  been  separated  out. 

In  equation  (19);  and  the  quantity  |i'  nov  r^resents  a bulk  viscosity 
which  pres\unably  takes  account  of  all  Intra -molecular  processes.  In- 
cluding molecular  relaxation  idienomena  plus  any  "inherent"  bulk  vis- 
cosity (if  there  Is  any  such  quantity). 

Both  theory^  and  recent  experiment^  support  the  assuiqptlon 
that  (1*  ■ V*  V 0 for  argon.  The  attenuation  Is  definitely  clas- 
sical, and  therefor*  one  would  expect  the  e creaming  to  dq;>end  only  on 
the  classical  mschanlsms  of  shear  viscosity  and  heat  conductivity. 

The  constants  for  argon  are  very  well  known,  and  substitution  Into 
equation  (19)  gives 


b - 2.25 
argon 


The  experimental,  streaming,  aetermlnatlon  of  b*  for  argon  was 
performed  by  observing  the  axial  values  and  vi202  from 


3U 

the  source.  Equations  (2^)  and  (26)  of  the  next  section  are  used  to 
deteroflne  the  geosetrlcal  o«istant  of  pr<^rtlonallt7  betveen  1*2.02 
and  UgQ^y  vhich  le 

U202  - .00940  KPiog  r/  • .00940  (20) 

P ^ cj 

o o 

The  critical  Irfornoatlon  needed  for  the  evaluation  of  h*  for 
argon  Includes  the  calibration  of  the  A12P  microphone  In  argon.  Flnt 
a WB  64QAA  microphone  vas  calibrated  by  free-field  reciprocity  at 
183  EC  in  air  to  give 

« - 97.9  db  re  1 volt  per  dyne  / cm^. 
nien  a flat  (±  1 db)  region  of  the  field  of  the  BaTlO^  source  vas 
located  and  the  64QAA  and  the  traxxsverse  XDP  alternately  placed  In 
this  region  to  give  the  relative  response.  Taking  account  of  the 
dlfferraices  betveen  the  gains  of  the  two  pre-anpllflers.  It  was  found 
that 

Hus  ■ %,0  - 17.9  db. 

Finally,  since  the  end  of  the  ADP  was  approximately  one -half  wave- 
length In  thickness,  It  was  necessary  to  observe  the  dlfi’erence  in 
sensitivity  of  ASP  In  air  vs.  in  argon.  (The  wave  length  In  air  Is 
approximately  10^  greater  than  In  argon.  The  beam  diffraction  was  ex- 
perimentally found  to  be  Independent  of  this  small  difference  In  wave- 
length which  was  as  expected.  However,  the  relative  else  of  microphone 
to  wavelength  made  the  microphone  very  sensitive  to  small  changes  In 
wavelength.)  With  the  640AA  at  the  end  of  the  tube  to  monitor  the 
field  Intensity,  change  of  v^ivm.  from  air  to  argon  was  accompanied  by 
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a decrease  In  AIP  sensitivity  of  4.2  dt.  The  resulting  sensitivity 
of  the  ASP  In  argon  vas 


argon  * db  re  1 volt  per  dyne/cn^. 


With  the  ajtollfloatlon  of  52.5  db  and  the  Mcrophone  Com- 
plement plxis  Pre •amplifier  gain  of  35.8  db.  cm  output  of  1 volt  na 
meant  an  AUP  oul^t  of  >88.3  db  re  1 volt.  This  represents  an  acous- 
tic field  of  31.7  db  above  1 dyne/cn?,  or  an  acoustic  rms  pressure  of 
38.5  dyne/co^.  The  aocon^anylng  streaming  reduced  to  one  volt  oul^t 
vas  0.0268  cm/sec.  These  values  and  the  constants  for  argon  give 


2 c 5 

Po  % 


0.009406^r  ®(P?o2/2 


2.34 


idilch  differs  from  the  classical  or  attenuation  value  of  b by  only 


At  this  Juncture  it  is  vortlndille  to  restate  the  assuoptions 
idilch  were  made  In  the  theoretical  solution,  and  which  seem  to  be 
supported  by  this  e3Q>erlmental  check.  Ihese  assusptlons  were: 

1)  that  the  Wavier *Stokes  equation,  the  equation  of  con- 

2 

tlnulty  and  the  acoustic  equation  of  otate  p^  ■ ‘^(aPi 

provide  a valid  basis  for  second -order  studies  of  streaming; 

2)  that  thermal  conduction  effects  may  be  Introduced  Into  the 
Sckart  streaming  equations  in  the  same  manner  as  they  appeeir 
in  attenuation  formulas; 

and,  mor*e  particularly  in  this  experiment, 


This  is  the  average  of  four  separate  series  of  observations  (a 
total  of  twenty  timings)  using  both  lycopodium  and  tobacco  smoke  as 
tracing  particles,  in  two  diffezvnt  gas  saoiples.  The  maxlraim  devia- 
tion of  the  four  values  from  the  mean  value  presented  above  was  6f». 


36 

3)  that^  was  the  oal7  effectlra  aouroe  of  vortlclty; 

4)  that  the  dtrergenoe,  attesoatlon  and  departure  from  plane - 
wave  character  of  thla  aound  bean  penal tted  the  aaauaftlon 

-1  " aln  (iDC  - ut) 

rather  than  the  nore  accurate  (and  Tlrtually  unaanageable) 

deacrlptlon  - ax  - Jut 

Pi  ■ Pi^**#*>e 

The  validity  of  aasuaptlcn  4}  la  dlacuased  In  Appendix  B. 

4)  Attenuation  of  Axial  Strean  Velocity  (Quarta  Wind) 

In  the  earlier  explanations  of  the  "quartz  vlnd,”  the  fluid  flov 
vas  ascribed  to  an  111  •defined,  local  reotlfyins  action  due  to  the 
vibrating  crystal.  It  vlll  be  shovn  here  that  Sokart'a  theory,  vhen 
properly  applied  to  the  actual  acoustic  field,  gives  a quantitative 
explanation  of  the  observed  decrease  in  axial  streaming  velocity  as  one 
recedes  from  the  soiuoe. 

The  variation  In  axial  stream  velocity  vas  observed  at  various 
distances  from  20  cm  to  72  cm  from  the  source.  Three  runs  at  differ- 
ent fields  Intensities  vere  made,  and  the  linear  relation  between  In* 
tensity  and  stream  velocity  (Figure  vas  used  to  refer  all  valties  to 
the  same  acoustic  field  Intensity.  By  defining  the  streaming  speed  at 
z ■ A In  db  as 

10  loSlO 


and  using  the  reference  velocity  U2op  * 1*00  cm/sec,  the  axial  stiream- 
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log  speeds  at  different  dlstanoes,  s,  from  the  source  oan  ha  plotted 
In  dh  re  1. 00  cm/sec.  This  Is  shovn  In  Figure  11.  The  change  In 
speed  per  unit  distance, 

s . U/z)  10  logjp  f.^')  (aa) 

Is  defined  as  the  axial  streanlng  attenuation.  The  value  of  S Is 
observed  from  Figure  11  to  be  0.16^  db/cm. 

In  order  to  determlxM  the  causer,  of  this  stream  attenuation,  the 
equation  tdilch  predicts  the  axial  streaming  velocity  may  bo  vrltten, 
for  a given  axial  position  A,  as 


2QA 

or,  for  position  B,  as 


®A^10A^  C 


(23) 


'^iOB 


The  letter  G represents  a factor  vhlch  takes  Into  account  the  variation 

of  UgQ  vlth  geoBBstrloal  form  of  the  acoustic  field,  at  constant  axial 
2 

Intensity  . C Is  a constant. 

Forming  the  ratio,  taking  the  log  to  the  base  ten,  and  dividing 
by  the  separation  c « there  results 

U/z)  10  lx«ja  ^ . U/Z)  10  logjo^V^)  ♦ U/2>  20  log^  ^ 

2 *8  ♦ a (24) 

Equatioa  (24)  describes  the  axial  streaming  attenuation  In  db 
per  unit  distance  as  caused  by  an  equivalent  attenuation  due  to  a 
divergence  of  the  acoustic  beam,  6 , and  a real  acoustic  Intensity 
attenuation,  a , due  to  absorption  of  acoustic  energy*  Such  a pre- 
sentation relies  on  the  validly  of  using  a psrturb;^tlon  of  the  Bokart 
solution  to  take  account  of  the  changing  field  at  different  distances 
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from  the  source.  This  point  is  disoussed  In  Appendix  B. 

To  calculate  the  term  due  to  ohange  of  the  acoustic  geometry, 

the  acoustic  field  vas  recorded  at  26  om,  cm,  and  12  cm  from  the 

source.  At  each  distance  from  the  source  the  microphone  probe  vas 

moved  through  four  different  diameters  and  these  responses  averaged  to 

give  the  average  distributions  shovn  In  Figure  9.  All  records  vere 

taken  with  the  recorder  potentiometer  adjusted  so  that  the  recorder 

deflected  full  scale  for  the  pressure  on  the  axis,  But  these 

axial  pressures  woul^l  not  be  equal,  even  on  the  assuiqptlon  of  no  at> 

tenuation.  Bather,  the  energy  traxismitted  through  any  cross  •eectlon 

of  the  tube  should  be  the  same  if  there  is  no  attenuation.  So 
r**o  2 

1 rPj^Q  dr  vas  cospited  graphically,  from  curves  such  as  In  Flguire 
7,  for  these  three  different  cross -eections  to  give: 

)o  rf-io®  . 0.0268 

- O-Ol"*? 

- 0.0073 

To  "nomallze"  vith  respect  to  the  intensity  at  12  cnu 


1,80 


(25) 


Bsxt,  the  axial  stream  reloclty  v&s  calculated  by  gmphlcal  In- 
tegratlon  of  ( Q P^(r)  dr  from  curves  such  as  In  Figure  8,  with 


the  results: 


uo 

T1202  ■ 0-0345  - 2.83  db  re 

ugo5  - 0.0261  r/  • 1.62  db  re  Ug^^  ^ggj 

ugo7  - 0.0180  E»io7^  **o^  - 0 db  re  Ugo^ 

Oliese  velocities  are  plotted  In  Figure  12.  The  attenuation  of  the 
axial  stream  velocl^  duo  to  divergence  of  the  acoustic  beam  Is  ob« 
sexved  to  be 

8 it  0.062  db/cm. 

The  term  a ■ (l/z)  20  1o6,q  is  the  sound  Intensity  attenu> 

*^1CA 

atlon  per  unit  distance.  The  value  a « O.078  db/cm  for  air  at 
185  KC  given  by  the  classical  theox7  Is  generally  conceded  to  be  too 
low,  because  only  shear  viscosity  aixd  thermal  conductivity  are  con* 
sldered  as  soxzrces  of  energy  loss.  Beyond  this,  there  Is  general  and 
great  disagreement  about  the  magnitude  of  the  correction  to  the  clas* 
slcal  formula.  The  vo3dc  of  Slvlan^^  Is  frequently  quoted  In  the  lit- 
erature. He  gives  a « O.llS  db/cm  for  dry  air  axid  a ■ 0.144  db/cm 
for  air  vlth  l,26f>  H2O  molecules,  but  states  that  his  data  are  "subject 
to  an  r.m.s.  spread  of  the  order  of  7;5^.”  More  recent  worh^^  leads  to 
a slightly  lover  value  for  moist  air.  The  e:Q)erlments  described  In 
this  secti<Ki  were  performed  using  air  of  approximately  1^  BgO  mole- 
cules. The  val\ie  of  the  acoustic  Intensity  attenuation,  a , idien 
added  to  the  divergence  attenuation  factor,  8 , sho^ild  lead  to 
the  axial  streaming  attenuation,  Z • 


I 

f 


1^  - 


f 

i 


Source 

i HgO 

a 

db/cm. 

6 

db/cm. 

Z - 

0+5 

db/cm. 
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Dlscnp 

Slvian  Moist  Air 

1.26 

0.144 

0.062 

0.208 

+ 25 

Penn  State  Moist  Air 

1.00 

0.126 

0.062 

0.188 

+ 14 

Slvian  Dry  Air 

0 

0.118 

0.062 

0.180 

+ 9 

Classical  Dry  Air 

0 

0.078 

0.062 

0.140 

- 15 

nie  values  In  the  fourth  coltum  are  coeqpared  vlth  the  observed  axial 
streaming  attenuation,  » O.165  db/cm. 

Susmarltlng:  the  observed  axial  stream  attenuation  of  0.16^ 

db/em  in  moist  air  can  be  accounted  for  by  the  change  In  the  acoustic 
field  Intensity  and  geometry  as  the  sound  wave  proceeds  avay  from  the 
source.  The  effect  Is  not  a local  one  at  the  source,  but  takes  place 
in  the  medium. 

C.  DSPSHDEITCE  OF  STSEAMIK&  OH  TBS  MS3>I0M 


Having  verified  that  Introduction  of  thermal  conductivity  Into 
the  definition  of  b*  gives  an  excellent  agreement  between  b*  obtained 
by  streaming  and  the  value  from  attenuation  measurements  for  argon, 
one  looks  for  a similar  situation  for  other  gases,  too. 

Continuing  the  analysis  of  b*  as  Including  all  attenuatlve 
mechanisms,  the  non  •classical  part  of  equation  (19)  may  be  written 


NC 


2 pea  jjQ 

k2“ 


(27) 


vher9  the  classical  effects  of  heat  conductivity  and  shear  viscosity 
have  been  subtracted  out. 


T^:e  mechanism  of  molecular  relaxation  has  been  sucoessflil  in 
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explaining  the  non>olu8lcal  absorption  for  aany  gases,  and  the  aqpll- 


tude  attenuation  In  nepers  per  unit  length  by  this  prooess  Is 

2 


M 


Y -1  4x0 


°oY 


(28) 


^10 


4-  U 


^ere  b Internal  specific  heat 

Cy  • specific  heat  at  constant  toIuisb 
Oq  * speed  of  sound 
Y B ratio  of  specific  heats 
kxo  * 2nfxo*  angular  relaxation  frequency 

u • 2irf  ■ angular  acoustic  frequency  at  vhlch 
molecular  attenuation  Is 

The  process  described  by  equation  (28)  considers  that  there  Is  a 
relaxation  frequency  for  a glren  lntra>molecular  process  (rlbratlon  or 
rotation)  such  that  idien  the  fluid  Is  excited  by  sound  ware  at  this 
frequency,  thex«  Is  a maximum  acoustic  attenuation  per  varelength. 

This  energy  loss  Is  due  to  the  mcleoules  being  set  Into  vibration  (or 
rotation)  and,  because  of  the  time  lag  necessary  for  taking  up  and  re- 
turning this  ensxgy,  the  energy  comes  back  to  the  translational  mode 
out  of  phase  with  the  sound  vave,  and  Is  thereby  degenerated  Into 
heat.  ?or  lover  or  hl^wr  acoustic  frequencies  than  this  relaxation 
frequency,  the  loes  of  enezgy  per  cycle  Is  less.  Presumably  such  a 
process  vould  shov  up  as  the  X or  p*  of  the  oaorosooplo  hydrodynaolo 
equations  (l)  or  (2),  From  the  results  of  section  II-B-3  It  Is  clear 
that  the  monatomic  gas  argon  has  zero  bulk  vlncosity.  This  section 
shall  study  the  question  of  vhether  there  Is  an  Inherent  bulk  viscos- 
ity In  the  gases  dry  nitrogen  and  moist  air,  or  irtwther  the  bulk 


kk 


\ 

i 


viscosity  l8  COH|pO80d  t^d  CCuwX*iwUlilO&  i.nt2?&* 

BDlecular  processes. 


S^uatlcms  (27)  axid  (26)  lead  to 


3o'*M 


(Y-1) 


kio^  + 


where  p ■ ataoeph8x*lo  pressure 

(1*  ■ hulk  viscosity  due  to  molecular 
M relaxation 


The  value  of  » Is  seen  to  reduce  to 


Cl  Y-1 


■ a constant  for  u < < k-t 


• ®1  (Y-i)  ^10 

»*•  ■‘h-'oT 


for  w > > k. 


So  the  value  of  ii'jj  Is  zero  for  very  high  freq^uencles  with  re- 
spect to  the  relaxation  frequency  and  climbs  through  the  relaxation 
frequency  to  reach  a steady  value  for  all  frequencies  much  less  than 
the  relaxation  frequency.  Schematically: 


Nitrogen  Xj^q  * 

The  streaming  experiment  was  conducted  at  l8^  KC  and  so  one 
would  expect  contributions  to  from  both  vibrational  relaxation 
which  occurs^  strongly  In  the  region  15  KC  to  25  KC,  and  from  rota- 
tional relaxation  which  occurs  In  the  b50  megacycle  region.  Schmldt- 
mtiller  observed  attenuations  from  KC  to  100  KC  In  nitrogen.  He 
found  that  the  exoese  (above  classical)  absoiption  could  be  explained 
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equalljr  veil  as  being  due  to  either: 

a)  a vibrational  relaxation  frequency  at  25  KC,  or 
B)  a vibrational  relaxation  f^quenoy  at  15  ICC,  supcrlqposed  on 
an  "excess  classical"  absorption  of  15^. 

On  the  basis  of  his  hypothesis  A),  and  using  his  constants: 


there  results 


- 0.0029 

• 5.0 


Y - 1.40 
p °1  (Y“1)  ^10 

^10  * ^ 


ho  - 25  -C 

f ■ 185  KC 


0.151 


t^ere  Is  the  bulh  viscosity  contribution  due  to  nolecular  vlbra* 

tlon. 

Or  using  his  hypothesis  B),  In  vhlch  the  vibrational  relaxation 
frequency,  f^^.  Is  assuiasd  to  be  15  EC,  there  results 


MT/ 


0.10  and  li'  Excess  classlcal/p  ■ 0.26 


Talues  of  • 


from  attenuation  or  streaming 


experiments  are  obtained  by  subtracting  b . • 4/3  * ■ 1.85 

Class 

from  the  value  of  b*  as  defined  by  equation  (l9).  These  values  are 
all  for  the  frequency  l85  EC.  They  are  tabulated  belov; 


Table  I 


b* 

Source 

2.00 

0.15 

Schmldtmuller, 

hypothesis  a) 

2.18 

0.33 

Strszialng 

2.21 

0.36 

S chmldtmuller , 

hypothesis  B) 
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Clearly,  the  streaoicg  experiment  support  hypothesis  B)  In  prefer- 
ence to  h^othesls  A).  But;  at  the  sane  time,  the  1%  correction  to 
classical  attenuatl<»i  nov  demands  “une  raison  d'etx«.'* 

One  vould  expect  that  when  Schmldtnuller  adds  1^  to  his  clas- 
sical attenuation  he  Is  taking  account  of  the  contrlbutlc»a  of  rota- 
tional relaxation,  so  that  classical  labeled 

H 4-,  Today,  following  the  work  of  Zmuda-^'  and  Parker,"^  this  cor- 
rectic’;*  should  be  ou  a better  experiiflental  basis . The  Parker  at- 
tenuation ejgpexlment  used  a progressive  pulse  method  with  frequencies 
from  60  EC  to  70  EC  In  the  range  ^ to  30  megacycles  per  atmosphere. 

The  work  by  Zigj.da  was  perfonaed  by  Interferometer  at  the  frequency 
3 megacycles  per  second  in  the  range  3 to  109  megacycles  per  atmos- 
phere, and  shovs  the  beginnings  of  rotational  relaxation  towards  the 
higher  f/p  values.  Although  the  frequency  ranges  of  these  two  e:q^rl- 
msnts  overlap,  the  results  show  a discrepancy,  the  Parker  wuxk 
giving  the  ratio  of  experimental  at  iatlon  to  classical  attenuation 
as  1.2^,  idille  Zouda  gives  1.40.  The  lower  frequencies  of  both  ex- 
periments are  svifflciently  small  with  respect  to  the  rotational  re- 
laxation frequency  for  to  be  given  by  the  constant  of  equa- 
tion (30).  The  Parker  results  are  probably  the  more  reliable  of  the 
two  because  a)  the  validity  of  the  measurements  was  first  verified  by 
a successful  determination  of  the  attenuation  in  argon,  and  b)  the 
value  of  the  ratio,  experimental  attenuation  divided  by  classical  at- 
tenuation, Is  close  to  1.2  irtiioh  Is  the  result  to  be  e3q>ected  from 
classical  kinetic  theory  for  rotational  relaxation.  The  Patker  value 
of  the  rotational  relaxation  frequency  will  be  used  to  calculate  the 
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value  of  (A'}/|p. 

n'  MR/ti  = -2  ^ = 0.31 

li  Cv  kio 


where  = 2.88  x 10^  eec"^ 

= B 

= 5/2R 
Y-1  = O.UO 

p = 1.013  X 10^  dyne/cn^ 


This  value  la  equivalent  to  175b  of  = I.85  and,  therefore, 

reinforces  the  argument  for  hypothesis  B). 


One  may  conclude  that  the  results  of  Schmid tculler,  Parker  and 
this  streaming  experiment  concur,  and  that  the  attenuation  of  soiind  in 
dry  nitrogen  is  due  to  the  classical  lass  mechanisms  of  heat  conductive 
Ity  and  shear  viscosity,  plus  losses  due  to  a vibrational  relaxation 
frequency  at  approximately  I3  EC  and  a rotational  relaxation  frequency 
at  approximately  4^0  megacycles.  This  streaming  experiment  provides 
no  evidence  of  an  Inherent  bulk  viscosity  which  is  distinct  from  relaxa- 
tion contributions. 


Moist  Air 

In  the  case  of  moist  air  of  1*3^  EgO  molecules,  there  is  avail- 
able the  data  of  attenuation  eiqperiments  at  Pennsylvania  State  College^^ 
and  the  correction  for  molecular  absorption  due  to  H2O  given  by  Xheser.^9 
The  value  of  b*  obtained  from  streasdag  falls  between  these  two  V8^.ues, 
again  supporting  the  argument  that  the  streasiing  experiment  is  another 
method  for  obtaining  attenuation  constants.  The  three  values  of  b^ 
are  given  below,  in  Table  II. 


III.  CONCLUSIONS 


The  vork  of  Sckart,  vhlch  conalders  the  vorticitjr  generated  in 
a vleeouB  medium  during  paaeage  of  a aound  beam,  haa  been  extended 
to  conalder  "atirface  aourcea"  of  vortlclty  vhlch  appear  vhen  the  aound 
beau  makee  contact  vlth  a a olid,  reatralnixig  aurface.  Aa  in  the  Eck- 
arc  'oork,  the  dynamic  ahear  vlacoaity  coefficient,  ^,  ia  aaauued  to 
be  Independent  of  the  denalty  changee  which  accompany  a aound  wave. 

However,  the  effect  of  a deneity  dependent  bulk  vlacoaity  coefficient 
la  inveatlgated,  and  it  ia  foruid  that  the  aourcea  of  vortlcity  are 
unaffected . 

A guided,  progreealve  wave  in  which  tube  radiua  < acoustic 

haa  been  atudied.  The  predominant  aourcea  of  vortlcity  in  this  case 
are  the  "surface  sources,"  and  the  resultant  acoustic  streaming  is 
controlled  by  these  "surface  sources." 

A radlally-sysinetrlc  progressive  acoustic  beam  haa  been  set  up 
in  a closed,  gas-filled  tube  to  cause  streams  driven  by  the  Eckart, 
volume,  source  of  vortlcity.  A careful  determination  of  the  acoustic 
pressure  at  all  points  of  a cross-section  of  the  beam  haa  been  used 
to  auceesatully  predict  the  stream  velocities  at  this  cross-section, 
for  argon.  The  work  vlth  argon  and  the  rosulta  of  further  experi- 
ments vlth  dry  nitrogen  and  Boiat  air  indicate  that: 

1)  the  Navler-Stokea  equation,  the  equaticn  cf  continuity,  and 

g 

the  "acoustic  equation  of  state," 

provide  a valid  haale  for  second-order  etudlee  of  etreamlng; 
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2)  theraal  conduction  effects  nay  be  Introduced  Into  the 
Sckart  etreanlng  equatione  in  the  sane  manner  aa  they  ap» 
pear  in  the  attenuation  formulas; 

5)  the  macroscopic  quantity bulk  viscosity,  is  made  up  of 
contributions  from  the  intra-molecular  processes  of  vibra- 
tional and  rotational  rej^txation. 
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Appefndlir  A 


STREAMIIIG  CAUSED  BY  A PLANE  PBOGBESSIVB  KAVE  IN  A CYIINDRICAL  OUBE 

The  purpoee  of  this  appendix  la  to  coopare  the  relative  magnl* 
tvdee  of  the  surface  aikd  volume  aourcea  of  streaming  for  a guided  plane 
progressive  vave  In  a cylindrical  tube. 

In  order  to  assess  the  relative  iaportanee  of  the  sources  of 
vortlclty  represented  by  the  three  terms  cc  the  right  of  equation  (^), 
it  is  convenleat  to  present  the  equation,  in  so  far  as  possible,  In 
tenns  of  the  first  order  acoustic  particle  velocity,  u^^,  or  its  curl, 

Using  (2)  and  (3),  gives, 

72^2  = -(b/co2)(a;Jj^/dtx77.Ui)  - [bv^/Co^]  [ (7lBi)  x77.Ui]  (32) 

-(1/Vq)7x(ui3Bi)  + (l/vQ)(dv/dp)o7X[pi(7xBi)] 


In  equation  (32)  the  first  two  sources  of  vorticity  have  cmbb 
from  the  tern  ^ . Identifying  the  second  as  "A"  and  cooparing 
with  the  third  term,  Sp,  in  equation  (32),  "A"  can  be  shown  to  be 
generally  negligible . 


Expanding:  A = [-bv^/cg^]  ((7:®j)  x72ujl] 

while  ^ f 

Therefore:  !^|/|a|  = (l/vo)/(vobk//c/)  . (*»A>)(>‘acous.Avlac)*^ 

where  . = 2ff  (2Vo/q)^/2 

Xacous  /*o  Xvlsc 
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Since  it  it  g«Mrally  tme  that  the  vieeoua  wave  length  la  aueh  leea 
them  the  acoustic  vave  length  In  fluids,  the  "A"  teiis  viU  henceforth 
he  dropped  and  vhat  ia  left  of  ^ vill  he  called  8^.  Then  the  vor> 
ticity  equation  nay  he  vritten: 


^^2=  +(-h/c^’^)(aUj^/5txv7*Uj^)  + (•1Ao)[7»(Ui^3^j^)  3 + (lA^,)(dv/dp)^7x[Pj^(7«j^)] 


3. 


Sr 


Sip 


(35) 


To  proceed  further.  It  is  neceaaary  to  specify  the  boundary  con- 
ditions of  the  particular  proiblem^  so  that  expressions  for  U]^  and 
nay  he  presented.  A search  of  the  literature  reveals  no  exact  solution 
to  the  Bavler-Stokea  equation  In  cylindrical  coordinates,  although  pre- 
llainaxy  equations  have  been  set  up  by  Eirchoff  and  reproduced  by  Ray- 
lel^^  for  the  ease  of  a guided  plane  progressive  vave  in  a cylindrical 
tube.  One  of  tbs  constants  of  integration  Is  evaluated  by  the  condi- 
tion that  the  particle  velocity  be  sero  on  the  tube  vails.  The  solu- 

1 as 

tlon  based  on  > k then  becosns: 

Axial  Velocity,  ^ A { J^[(-l+l)3r]  - [ (-l+i)8B  ] } (34) 

Radial  Velocity,  - [Ak(M)/(gp)3  {-  Jil(-l+l)Sr  ] + (r/R)Ji[(-l+i)0Bj} 

(55) 

vhsret 

S3  = frequency  of  the  acoustic  vave 

k » > aeouetlc  vave  nunher 

ftOOUUa  e 

R s tube  radius 


*?or  axasgple,  for  a nonatoedc  gae  ^^acous  **  **  **** 
period  of  the  sound  vave  is  long  cooffairsd  to  the  tins  betveen  atosdc 
collisions. 
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r e distance  frcrn  tuba  axis 
h ::  B-r  > distance  from  ttd>e  vail 

p = ^(0/2  V 

V = ti’/ps  kinematic  shear  viscosity. 

An  approximate  form  of  this  solution,  vhlch  holds  for  r > 0 and 
for  0r>l,  for  a vave  propagated  In  the  -fZ,  1^,  direction  and  of  vel* 
oclty  coB(kZ-ut)lz  on  the  axis  is: 

Wi  = ^ cos(kZ-«t)-e“^^  coB(kZ><iit>0h)} 

Ui  = [(Uj^Qk)/>/20]  { (r/fi)coe(kZ-wt4-w/4)-e’^^  cos(kZ-ut-0h-ftr/U)) 

The  assumed  solution  is  seen  to  have  a radial  ss  veil  as  an  axial  con- 
ponent,  both  vith  a precipitous  drop  in  the  particle  velocity  through 
the  thin  boundary  layer  near  the  vail.  0 is  a large  number  (e.  g.  for 
f » 1000  cps  in  air,  n = 2 x 10"^  poises,  3-200,  cm"^).  Further,  the 
radial  ccspcnent  is  smaller  than  the  axial  by  virtue  of  the  factor  k/3 
vhlch,  for  the  exsople  given,  vould  be  about  10~^. 

Expanding  all  vcrtlclty  sources  in  cylindrical  coordinates  r,  9, 

Z vhsn  there  is  axial  syaastry  all  derivatives  vith  respect  to  e 
may  be  set  equal  to  zero.  Then,  denoting  the  unit  vector  in  the  6 
direjtloo  by 

8e  » -(b/c/)(9Ji/atx7V-‘Si)  - i -UtLOrz  + UA)^z+  "zz? 

(37) 

h • -a/v^)vx(;^«^) « igC-i/v^)  { -w 

+ Uj.Wr> 
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Sj  - = -(i2/po){-SPi/ar  [U„-W„4Uj/r-Wj./r3  -3Pi/«a 

[Ujj-Wjj.]  + Pj^  L~^zzz'^2zr~^^!Tz^rrr*^rzA‘'*'  ^TzA^'*^rr/**“'^r/*^  ^ ^ 


vharrj  the  ahoi’t-hand  notation  Uj.^  = 3®U^r8B  , Wj^  li  3%^r^  etc. 

Using  equations  (36)  for  the  evaluation,  a large  nus2>er  of  quad* 
ratio  terns  are  fcmed.  Of  these,  only  the  ones  with  a non-zero  tine 
leverage  val\ie  are  of  Interest  In  a discussion  of  streanlng;  the  second 
hanacmlc  tarns  are  dropped. 


For  the  conditions  3>l/r>k  there  results: 

Sg  - -l2(t0’k2u|Q/co){-e“^^  + e"^^[  slaph  + cos^h]  + 

Sg  » -l^Cpkufo/v^)  [ (r/fe)<sln3h+cos3h)]  } (58) 


^ . -IgCpkuioAo) 


e-23h  - 1 e'^^(sln3h-coa3h)\ 


{^rr  01- 

It  Is  noted  that  the  terminology  adapted  for  the  sources  of 
vocrtlclty  Is  supported  by  the  presence  of  h (vhlch  Involves  hulk  vis- 
cosity) only  In  the  volume  term,  The  sources  of  vortlclty  are 
similar  In  their  variation  vlthln  the  housdazy  layer  hut  only  the 
volume  source  has  a term  vhlch  takes  on  significant  values  even  avay 
from  the  vail.  It  Is  of  Interest  to  coopars  the  orders  of  magnitude 
of  th>:  sources  of  vortlclty  and  It  appears  that^ 


I Surface  Sources 


f£_  - iAwt) 
h(ovo 


jVolums  Source  | 

Where  t « (4/5)vo/C(,^  « "viscous  relaxation  tins"  as  used  hy  Eckart,  and 
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la  of  the  aaae  order  of  nagnltude  aa  the  tin  between  coUlaloia  for  a 
Bomt(»slc  gas. 

Since  the  nagnltude  of  the  voluae  aource  of  vortlclty  la  propor- 
tional to  frequency  and  the  surface  tern  are  not,  one  locdca  for  a 
frequency  at  vhlch  the  volume  aource  would  become  of  Inportance  for  a 
sound  beam  which  fills  the  tube.  Apparently  this  would  occur  when  the 
frequency  approaches  (1/t).  However,  it  should  be  pointed  out  that 
the  equations  (36)  un  which  this  conclusion  la  based  are  valid  only 

lf\.  . But  within  these  limits,  the  surface  sources  of 

vise  acoun  ' 

vortlclty  appear  more  Important  than  the  volume  source.. 

But  the  answer  to  the  obvious  question  of  whether  the  streaming, 
say  on  the  axis.  Is  controlled  by  the  surface  sources  and  In  the  same 
ratio  as  above,  req\ilres  an  Integration  of  equations  (38) , A vector 
potential  defined  by  V^A2  > -B2  (vbere  A2  Is  a solenoldal  vector) 
la  proposed.  Then  s 73^2  for  an  incmm^resslble  fluid,  (Hote:  It 

la  assumed  that  IU2I  <<C^  so  that  the  streaming  flow  la  Incoo^resslble.) 
For  the  geometry  considered  here,  U2  « f(r)  Is  and  A2  » g(r)l2  where 
Ig  and  1^  are  the  unit  vectors  in  the  6 and  ^ (axial)  directions, 
respectively.  The  proposal  Aq  » r^  yields  the  solution  to  the 
homogeneous  equation,  A2  » (ar  + br^)l2.  It  Is  observed  tbat  the 
major  sources  of  vortlclty  are  near  the  wall  and  sc  the  assxmption  Is 
made  that  only  the  variation  of  A2  ^ ^ direction  near  the  wall 

Is  significant  In  determining  the  particular  Integral.  Therefore,  for 
the  uoterminatlcn  cf  the  particular  Integral  -»-V^(-7?A2)  Is  replaced 
by  d^Ae/dh^  a d^Ae/dr**,  After  integration,  the  two  constants  are 


evaluated  by  tbe  two  conditions; 
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(a)  W2  s 0 at  the  tube  wall 

(b)  W(r)  r A'  *0  for  mas  cootiaulty  through  a croaa  aection 
o 

of  the  tube. 

The  Poisaoc's  equation  la  aolved  flrat  with  aa  the  driving 


tem  to  yield: 


Wga  « + [bk^ufoA^^Co)]  {[k''0R^/5}tr‘*AB‘*  - r‘'/3H=^  t l/l2]  -l/2[l-2r‘'Al'^]} 


.2ots5  /e  “If  _ —2 


-2iu2i 


Then  the  Polaaon's  equation  la  aolved  vlth  8^^  ae  the  aouree  of  vor> 
ticlty  to  give: 

Wgj  » • ('^loA^o^  (l-2r^/R^)  ( 


Finally,  both  surface  aourcea  of  vortlcity  Sjj  + are  aet  up  aa  the 
driving  tema,  with  the  reault: 

%+2T  ‘ ' (3'»1oACo^  (l-2r^/R^) 

Only  the  etreamlng  outside  of  the  boundary  layer  la  reported  above. 

Study  of  these  three  reaulta  reveals  that  the  ratio  of  the  ng- 
nitxtdeo  of  the  two  types  of  atreaas  at  the  tube  eucls  is  approxlaately: 


I i 
•I  : i 


! . i 


'^2B42T  _ ®o 

w;: — 

2e  o 


and  so  the  atreaadng  velocity  in  this  case  la 


clearly  doalnated  by  the  surface  sources.  However,  this  ccscluslon 
Boat  not  be  extended  to  the  case  of  a sound  beam  (R>^an  ) which 


Since  kP  < 1 and  RS  > > 1,  It  la  possible  that  the  first  term 
In  the  Wgg  brackets  may  be  greater  or  less  than  the  second  term,  and 
could  even  result  In  W2Q  being  zero  on  the  axis . 


r 


6l 


does  not  fill  t&«  triibe.  Tbon,  oa  Bekart  points  out,  the  voluoe  source 
vould  generate  extr«Mly  great  stream  velocities  as  the  ratio  (beam 
radlus/tube  Badlus)'>0.  It  might  be  ejected  that  an  experimental 
sound  beam,  such  as  described  earlier  In  this  paper,  might  be  strong 
enou^  at  the  tube  valla  to  cause  a surface  source  of  vortlclty  vhlch 
vould  dominate  the  flow  bvit  this  Is  not  the  case.  Forming  the  ratio 
of: 


l“aoJ 


_0 
Bekart 


,OSb]^-  p2 


P p 2c  3 1 axis  ° 


to  |u_ 


20 1 
Surface 


«10*^ 


P 2 

*^1  vail 
4p  ^ 


there  results: 


Eck 

l'»20l 

Surf. 


itvi>2_ 


(15) 


(55) 


It  Is  clear  that  for  a real  sound  beam  In  vhlch  kr^  > > 1 and  > 

^1  vail'  the  Bekart  source  of  vortlclty  controls  the  axial  streamings 

Equations  59'’4l  present  another  fact;  that  S^,  In  general,  causes 
axial  flow*  avay  from  the  sound  source  vhlle  8p  and  cause  streamii^g 
avay  from  the  sound  source  at  the  vails,  returning  towards  the  source 
at  the  axle.  The  apparent  contradiction  Is  reconciled  vhen  cne  cen- 
alders  that  3^  Is  concerned  vlth  energy  loss  and  transfer  of  acoustic 
momentum  Into  streaming  momentum  vlthln  the  medliim.  The  surface  tersm, 
on  the  other  hajad,  take  account  of  Um>  transfer  of  momentum  as  the 
sound  energy  is  dissipated  vlthln  the  boundary  layer  at  the  tube  vails. 
A final  cbservatlon  Is  that  tbs  streaming  velocity  prediction  Is 
and  ^ are  used  as  driving  sources.  The  cause 


Increased  vhen  both 


of  tMB  inerease  is  recogalzed  If  the  tJiw  averago  of  8^  Is  written 


es: 


<h>  ■ pi 

o ® 

-> 

Then  shove  up  as  dependent  on  the  coe^resslhlllty  of  the  medium,  as 
well  as  the  Irrotatlonallty.  This  contrasts  with  Sj^  which  Is  a func- 
tion of  the  irrotatlonallty  of  the  medium,  only. 


APP2IKDIX  B 


VALEDITT  OF  TBS  ATTEKUaTIOII  APPROXIMATION 

The  Sckart  solution^  which  Is  the  baala  of  this  experlaental 
vcrh,  requires  the  Inclusion  cf  viscosity  In  the  force  equations. 

Later,  for  the  special  ease  of  the  "Ideal"  beam,  the  effects  of  bath 
viscosity  and  diffraction  on  the  flrst»order  acoustic  pressure  are 
coosidered  to  be  nesUglble,  and  the  equation  of  the  unattenuated, 
non-dlvergent  sound  beam  Is  written  as  « P(r)  sin  (kz  - a>t) . FuT' 
ther,  Eckart  has  assumed  that  there  are  no  end  effects.  The  result- 
ing solution  describes  the  streaming  as  Independent  of  longitudinal 
position,  z. 

Ezperlstentally,  since  the  streaming  at  the  center  of  the  face 
of  the  sinusoidally  vibrating  disk  must  be  zero,  there  Is  a region 
Mar  the  source  where  the  axial  stream  velocity  rises  rapidly  from 
zero  to  a maximum  value.  Slsdlarly,  the  axial  stream  velocity  must 
drop  to  a zero  value  at  tbs  fixed  end  plate.  Within  these  two  regions, 
the  streams  have  strong  radial  coe^onents  and  the  Eckart  solution  can- 
not be  expected  to  be  valid.  Telescopic  observation  shewed  that  the 
streaming  was  predominantly  longitudinal  in  the  region  from  20  cm  to 
7^  cm  from  the  source  In  this  eiq^riment,  and  all  data  were  taken 
within  this  region.  The  distance  from  source  to  absorbing  end  plug 
was  about  100  cm. 

Evan  this  central  region  does  not  fulfill  the  conditions  of  the 
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Xekart  Bolutlou  bee«u«e  of  attenuation,  dlverg«ice  and  noa>planar 
character  of  tbe  cound  vava.  The  latter  two  eondltlans  vould  he  dlf> 
floult  to  take  account  of,  and  are  probably  not  aa  Influential  aa  the 
attenuation  In  the  praaent  experinent  (aee  table  on  top  of  page  42 
vbere  the  nea-planar  effect  la  Included  In  the  geoasetrlcal  divergence 
term  6).  The  attenuation  could  be  deacrlbed  by  uaing  p^^  * Pj^(r)  aln 
(kz  - (ot)  e~^‘  vhere  a la  the  amplitude  atisnuntlon  coefficient . 
Subatltutlon  Into  ^ of  the  vortlclty  equation  (5)  leads  to 


v2b  _ ^*^4  d[Pj^2(y)j 

2 2p^2c^4  dr 


(42) 


a vector  Polaaon'a  iQuatlon  whose  aolutlon  vould  ahcw  that  the  atream> 
Ing  velocity  for  an  attenuated  bean  vlll  have  both  longitudinal  and 
radial  con^onente.  The  oagnltude  of  tbe  radial  coBi>onent  depends  on 
the  attenuation  a,  <•  beccHslng  zero  fcr  the  una.ttenuated  sound  bean. 

The  deterndnatlon  of  the  atreamlng  velocity  requires,  first, 
the  aolutlon  for  B2.  Solutions  to  the  honogeneoja  equation  are  aou^xt 
under  tbe  ecodltlona  of  ulal  symoetry,  regularity  at  the  axle  and 
dependence  on  both  r and  z.  Such  a solution  Is 

Rge  + Be"  (£r)  + Cr  +7*  (43) 

vhere  £ ia  & propagation  constant  for  vortlclty  and  74*  la  a poten- 
tial f;incticn  which  nay  be  Introduced  to  saiiafy  boundary  coaditiona. 
The  particular  integral  vould  not  be  aa  easy  to  find.  The  cooplete 

solution  for  U2  vould  require  the  integration  of  still  another,  more 

•*  •* 

ectiplex,  vector  Polaaon'a  Squation,  7<^A2  ■ vhere  A2  is  the 
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solezioldal  vector  potential  whose  cxurl  Is  the  streaaisg  velocity, 

Ug.  The  BoluticOL  viU  not  be  attenpted  here  but  It  should  be  nen> 
tloned  that  tbs  dctermlnatlosx  of  the  ccostants  of  Integration  and  the 
propagation  constant  for  vo&^iclty,  , would  have  to  be  evaluated 
frcm  a kncf^'ledgc  of  boundazy  conditions  on  velocity. 

Without  solving  this  problen,  scsaa  observations  nay,  neverthe> 
less,  be  ude.  It  is  noted  that  it  is  the  introduction  of  attenua> 
tlon  that  opens  the  way  to  solutlans  of  the  type  shown  in  equation 
(U3)  where  the  propagation  constant,  , appears  for  the  first  time. 
The  greeter  the  attenuation,  the  more  iiqportance  must  he  attached  to 
this  propagation  constant.  Vor  a hl^ily  attenuated  sound  beam,  the 
vartlclty,  generated  In  one  region  of  the  beam  can  be  es^ected  to 
propagate  Its  Influence  to  adjacent  sections.  Such  a propagation  of 
vartlclty  does  not  appear  in  the  Sckart  solxitlon  precisely  because 
the  acoustic  field  (and  the  reeulting  vortlclty  field)  which  he  treats 
has  no  dependence  on  z and  the  question  of  the  influence  of  one  cross - 
section  of  the  flew  on  the  adjacent  cross-sections  never  arises. 

In  lieu  of  an  analytic  eolutlon  to  this  prbblem,  the  author  has 
assumed  that,  within  the  limits  of  accuracy  of  this  experiment,  the 
streaming  in  each  section  of  the  tube  can  be  calculated  by  a knovlelge 
of  the  acoustic  field  at  this  same  cross-sectiem.  Physically,  this 
assumes  that  any  drag  or  drive  vhi  -n  nay  be  caused  by  adjacent  regions 
Is  negligible.  Otherwise  stated,  the  attenuation  is  ceneidered  to  be 
lerge  enough  to  affect  the  acoustic  field  but  not  large  exiough  to 
cause  measurable  ''coupling”  between  sections  of  the  streaming  field. 

As  one  proceeds  away  from  the  source  (In  the  middle  region  where  end 
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effects  are  negligible)  the  deoreaee  in  axial  etreaxdag  velocitlee 
ii  accounted  Tor  by  a relatively  eaall  **eeepage"  of  flow  into  radial 
cooponenta.  The  suba taut lat ice  cf  this  explanaticn  ie  entirely  ex- 
perisestal.  It  is  baaed  on  these  fantat 

a)  the  axial  and  off-axis  streaixlng  predictions  at  26  cm 
from  the  source  for  argon  and  nitrogen  coincide  (vithln  ^ 
vitb  predictions  based  on  the  observed  acoustic  field  at  26  cm 
frott  %hio  source* 

b)  The  axial  velocities  in  the  region  20  cm  to  73  cm  from 

the  source  are  predicted  from  a kncsrledge  of  (and  interpola- 
tion and  extrapolation  ^hroa)  the  observed  acoustic  fields  at 
26  cm,  ^0  cm  and  'J2  cm  from  the  source,  with  a naximum  error 
of  14^  (based  on  Pennsylvania  State  College  values,  loe . cit.). 
Vithin  the  limits  of  accuracy  of  this  experiment,  as  stated  above,  it 
is  felt  that  conclueloa  U)  of  page  36  is  valid. 
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